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INTRODUCTION 


A great  deal  of  work  has  been  accomplished  in  the  area  of  elastic 
laminated  effective  stiffness  or  microstructure  continuum  theories  and 
approximate  plate  and  beam  theories.  By  the  same  token,  little  has 
been  accomplished  with  viscoelastic  counterparts  to  these  theories. 

An  elastic  continuum  theory  which  included  effective  stiffnesses 
for  both  the  reinforcing  and  matrix  layers  of  a laminated  continuum 
was  developed  by  Sun,  Achenbach,  and  Herrmann  [1,  2],  The  continuum 
theory  was  utilized  by  Thomas  [3]  to  study  the  simple  thickness  modes 
for  laminated  media  with  layering  both  parallel  and  perpendicular  to 
the  plate  free  surfaces.  Sun  [4]  deduced  a two  dimensional  theory  for 
laminated  plates  from  the  three  dimensional  continuum  theory.  Veloc- 
ity correction  coefficients  were  introduced  into  the  two  dimensional 
theory  by  Thomas  [5]  and  flexural  and  extensional  vibrations  for  plate 

1C.  T.  SUN,  J.  D.  ACHENBACH  and  G.  HERRMANN  (1968)  Journal  of  Applied 
Mechanics , 35,  467.  Continuum  Theory  for  a Laminated  Medium. 

2J.  D.  ACHENBACH,  C.  T.  SUN  and  G.  HERRMANN  (1968)  Journal  cP  Appl ied 
Mechanics,  35,  689.  On  the  Vibrations  of  a Laminated  Body. 

3C.  R.  THOMAS  (1972)  Journal  of  Sound  and  Vibration,  23(3) , 341-361. 
Simple  Thickness  Modes  for  Laminated  Compos ffeTTaTe^Tal s . 

\.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  231-238.  Theory 
of  Laminated  Plates. 

^C.  R.  THOMAS  (1972)  Journal  of  Sound  and  Vibration,  25(3) , 407-431. 
Velocity  Corrected  Theory  of  Laminated  Plates  Applied  to  Free  Plate 
Strip  Vibrations. 
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strips  and  rectangular  pldtes  were  studied  by  Thomas  [6,  7]  according 
to  this  theory  and  compared  to  similar  results  from  effective  modulus 
plate  theories.  A microstructure  theory  for  an  elastic,  laminated 
composite  beam  was  developed  by  Sun  [8]  and  the  approach  utilized  in 
this  paper  will  be  followed  in  deriving  a viscoelastic,  laminated  com- 
posite beam  theory.  Thomas  [9]  showed  that  the  flexure  beam  theory  in 
reference  [8]  is  directly  obtainable  through  a simple  reduction  of  the 
existing  flexure  equations  for  composite  plates  [4,  5]. 

A continuum  theory  for  a viscoelastic  laminated  composite  was 
developed  by  Grot  and  Achenbach  [10]»  however  the  equations  developed 


^C.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  231-238.  Theory 
of  Laminated  Plates.  

5C.  R.  THOMAS  (1972)  Journal  of  Sound  and  Vibration,  25(3),  407-431. 
Velocity  Corrected  Theory  of  Laminated  Plates  AppTieB  to  Free  Plate 
Strip  Vibrations. 

6C.  R.  THOMAS  (1973)  Journal  of  Sound  and  Vibration,  31 (2) , 195-211. 
Extensional  Vibrations  of  Simply  Supported  Composite  Plate  Strips. 

7C.  R.  THOMAS  (1975)  Journal  of  the  Acoustical  Society  of  America, 
57(3) , 655-659.  Flexural  and  Extensional  Vibrations  of  Simply  Sup- 
ported Laminated  Rectangular  Plates. 

8C.  T.  SUN  (1971)  Journal  of  Appl ied  Mechanics,  38,  947-954.  Micro- 
structure Theory  for  a "Composite  Beam. 

9C.  R.  THOMAS  (1973)  Watery! iet  Arsenal  Technical  Report, 
R-WV-T-6-45-73.  Flexure  Equations  of  Motion”For~Larninated  Composite 
Beams . 

1 °R.  A.  GROT  and  J.  D.  ACHENBACH  (1970)  Acta  Mechanica,  9,  245-263. 
Linear  Isothermal  Theory  for  a Viscoelastic  Lamina  tel"  Composite. 


were  not  applied  to  any  problems  of  wave  propagation  or  vibration.  It 
is  certainly  theoretically  possible  to  start  with  the  equations  in 
reference  [10],  to  make  appropriate  series  expansions  and  derive  a 
plate  theory,  and  to  then  follow  reference  [9]  to  make  a direct  reduc- 
tion to  a viscoelastic  beam  theory.  However,  for  convenience  and  sim- 
plicity of  analysis,  the  approach  in  the  current  report  will  be  to  be- 
gin with  the  viscoelastic  Timoshenko  beam  equations  and  work  towards  a 
viscoelastic  laminated  beam  equation  in  the  manner  of  reference  [8]. 
With  somewhat  guarded  conclusions,  Stern,  Bedford,  and  Yew  [11]  have 
demonstrated  a definite  need  for  an  effective  stiffness  type  formula- 
tion for  viscoelastic  laminates. 

The  current  approach  to  obtaining  a viscoelastic  laminated  beam 
theory  will  be  a viscoelastic  development  which  mirrors  the  elastic 
development  given  by  Sun  [8].  Surprisingly,  the  real  difficulty  is  in 
obtaining  the  energies  for  a single  layer  modeled  as  a viscoelastic 
Timoshenko  beam.  The  most  pleasing  and  straightforward  development 
of  suitable  viscoelastic  Timoshenko  beams  results  from  a utilization 

®C.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  947-954.  Micro- 
structure Theory  for  a Composite  Beam. 

9C.  R.  THOMAS  (1973)  Watervl iet  Arsenal  Technical  Report, 
R-WV-T-6-45-73.  FI exure'Tquati ohs  oT  MoirTon  for  Laminated  Composite 
Beams. 

1 °R.  A.  GROT  and  J.  D.  ACHENBACH  (1970)  Acta  Mechanjca,  9,  245-263. 
Linear  Isothermal  Theory  for  a Viscoelastic  Laminated  Composite. 

1]M.  STERN,  A. BEDFORD,  and  C.  H.  YEW  (1971)  Journal  of  Applied  Mechan- 
ics , 38(2) , 448-454.  Wave  Propagation  in  Viscoelastic  Laminates. 
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of  viscoelastic  constitutive  relations  of  the  differential  form;  it  is 
these  equations  which  yield  a viscoelastic  development  which  closely 
mirrors  Sun's  [8]  elastic  derivation. 

THE  ENERGY  PRINCIPLE 

As  Sun  [8]  does  in  the  development  of  an  elastic  laminated  beam 
theory,  the  first  task  in  deriving  a viscoelastic  laminated  beam  theory 
is  to  formulate  energies  for  individual  viscoelastic  layers  in  terms  of 
the  Timoshenko  [12]  beam  theory.  In  the  past,  Lee  [13]  developed  vis- 
coelastic Timoshenko  beam  equations  for  viscoelastic  extensional  strain 
but  the  shear  strain  was  left  elastic.  Pan  [14]  extended  the  analysis 
to  include  viscoelastic  shear  strains.  The  current  objective  is  to 
develop  the  viscoelastic  Timoshenko  beam  equations  in  a form  more  suit- 
able to  the  development  of  a Viscoelastic  composite  beam  theory.  A first 
goal  will  be  the  development  of  a single  layer  energy  principle  suitable 
for  a direct  application  in  the  derivation  of  a multilayer  energy  prin- 
ciple. 

The  development  of  an  approximate  theory  such  as  for  laminated 
elastic  plates  has  originally  been  a two  step  procedure.  In  the  first 

^C.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  947-954 
Microstructure  Theory  for  a Composite  Beam. 

^2S.  P.  TIMOSHENKO  (1922 ; Philosophical  Magazine,  Ser.  6,  Vol . 43, 

125-131.  On  the  Transverse  Vibrations  of  Bars  of  Uniform  Cross-Section. 

nH.  C.  LEE  (1960)  Journal  of  Applied  Mechanics,  27,  551-556.  Forced 
lateral  Vibration  of  a Uniform  Cantilever  Beam  with  Internal  and 
External  Damping. 

1 4 

H.  PAN  (1966)  Jour.  Eng.  Mech.  Div.,  Proc.  Amer.  Soc.  Civil  Eng., 
213-234.  Vibration  ofa  Viscoelastic  Timoshenko  Beam.  ~ 

am*  " i i •m,nnnrra«lV~Mirfi 
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instance,  the  Mind!  in  plate  theory  [15]  in  its  first  order  approximation 
was  utilized  to  dewelop  a continuum  theory  for  laminated  composites. 
Tnen  to  obtain  a laminated  plate  theory  a first  order  approximation  is 
made  on  those  variables  in  the  continuum  theory  which  came  from  the 
zero  order  part  of  the  Mindlin  plate  theory  and  a zero  order  approxima- 
tion is  made  on  the  variaoles  which  came  from  the  first  order  part  of 
the  Mindlin  theory  as  in  Sun  [4]  and  Thomas  [5]  - this  explanation  will 
become  clear  shortly.  Now  in  developing  an  elastic  laminated  beam 
theory,  Sun  [8]  has  made  both  of  these  approximations  simultaneously 
to  obtain  a flexure  theory  for  laminated  beams.  Actually,  Thomas  [9] 
has  shown  that  the  flexure  beam  theory  is  directly  obtainable  from  the 
existing  flexure  plate  theory. 

The  current  objective  is  to  immediately  derive  a viscoelastic 
laminated  beam  theory  and  to  not  have  to  develop  a viscoelastic  lami- 
nated continuum  theory  first.  In  making  the  various  zero  and  first 


^C.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  231-238. 

Theory  of  Laminated  Plates. 

5C.  R.  THOMAS  (1972)  Journal  of  Sound  and  Vibration,  25(3) , 
407-431.  Velocity  Corrected  Theory  of  Laminated  'Plates  Applied 
to  Free  Plate  Strip  Vibrations. 

8C.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  947-954. 
Microstructure  Theory  for  a Composite  Beam. 

'C.  R.  THOMAS  (1973)  Watervliet  Arsenal  Technical  Report, 
R-WV-T-6-45-73.  Flexure  Equations  of  Motion  for  Laminated  Com- 
posite Beams. 

"1 8R.  d.  MINDLIN  (1955)  Signal  Corps  Engineering  Laboratories,  Fort 
Monmouth,  New  Jersey  '('Au-W'lTT)'.  An  Intro  due  t i on  "To  the  Mathe- 
maticnl  Theory  of  Vibrations  of  Elastic  Plates. 
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order  expansions  of  displacement,  terms  which  lead  to  an  extension 
theory  are  also  maintained  since  the  second  expansion  of  extensional 
displacements  leads  to  a flexure  term.  The  first  order  displa  ements 
which  will  result  in  the  Timoshenko  beam  equations  [12]  for  flexure 
as  well  as  an  extensional  equation  for  beams  are 
v(y,z,t)  = v(y,t)  - z<J>(y,t) 

w(y,z ,t)  = w(y,t)  - z<i>(y,t).  (1 ) 

the  zero  order  terms  in  (1)  are  v and  w and  a first  order  expansion  of 
these  two  displacements  results  in  the  expressions 
v(y,t)  = vk(y,t)  - zk  i|»a(y,t) 

w(y.t)  = wk(y,t)  - zk  H'a(y,t)  (2) 

where  the  subscript  a = 1,  2 will  later  denote  whether  a stiff  or  soft 
laminated  beam  layer  is  indicated  and  the  superscript  k which  layer 
pair  is  indicated.  While  absolutely  necessary  at  this  point,  the 
notation  in  (2)  jumps  into  the  laminate  notation  while  seeming  to  be 
at  the  single  layer  stage  of  development.  See  Sun,  Achenbach,  and 
Herrmann  [1]  or  Sun  [8]  if  clarification  is  required. 

Combining  equations  (1)  and  (2)  and  extracting  only  those  terms 
which  result  in  flexural  motion  results  in  the  displacement  relations 

]C.  T.  SUN,  J.  D.  ACHENBACH  and  G.  HERRMANN  (1968)  Journal  of  Applied 
Mechanics,  35,  467 . Continuum  Theory  for  a Laminated  Medium". 

8C.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  947-954.  Micro- 
structure Theory  for  a Composite  Beam. 

US.  P.  TIMOSHENKO  (1922)  Philosophical  Magazine,  Ser  6,  Vol  43,  125-131. 
On  the  Transverse  Vibrations  of  Bars  of  Uniform  Cross-Section. 
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v(.y,z,t)  = -z^a(/,t)  - z<j>(y.t) 
w(y,z,t)  = w£(y,t.) 


(3) 


where  v (y,t)  represents  the  gross  rotation  in  the  laminated  beam, 
wjj(y,t.)  represents  the  transverse  deflection,  and  c|>(y,t)  represents 
the  individual  layer  rotation.  The  various  displacements  and  rotations 
on  the  right  side  of  (1)  refer  to  the  individual  layers  which  will 
eventually  make  up  the  laminated  beam.  The  various  displacements  and 
rotations  on  the  right  side  of  (2)  represent  the  reduction  from  a lam- 
inated continuum  theory  to  a laminated  beam  theory;  thus,  from  conti- 
nuity of  displacement  and  rotation  at  laminate  interfaces,  it  is  clear 
that  the  notation  may  be  simplified  to  w(y,t)  = wj*(y,t)  and 
ijj(y,t)  - ^ (y ,t)  for  a = 1,  2 and  for  all  values  of  k.  Hence  with 
these  rotational  simplifications  in  mind,  the  final  form  of  the  first 
order  flexure  displacement  expansioi  is 

v(y,z,t)  = -zKjj(y,t)  - z<My»t) 

w(y,z,t)  = w(y,t)  (4) 

where  these  equations  are  valid  only  when  eventually  utilized  in 
developing  a laminated  beam  theory.  Equations  (4)  may  be  reduced  to 
those  for  a homogeneous  or  single  layered  beam  by  setting  ( y , t ) = 0; 
this  being  done,  equations  (4)  reduce  to  those  given  by  Brunelle  [16] 
for  flexure  of  a beam. 
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E.  J.  3RUNELLE  (1970)  J.  Composite  Materials,  4,  404-416.  The 
Statics  and  Dynamics  of  a Transversely  Isotropic  Timoshenko  Beam. 
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The  non-zero  strain-displacement  relations  are 
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The  non-zero  stress  equations  of  motion  which  pertain  to  the  problem 
are 


Jyz.y 


Jy>y 


= pw 


Jyz,z 


PV 


(6) 


From  the  appendix  and  equations  ( A- 1 7 ) the  constitutive  equations  for 
a special  case  of  the  standard  linear  model  are 


+ c ^)°yZ  = (2kG  + 2k  G ft)  Gyz 

(1  + C — )o  = (E  + E ~)  e 
3t  y 9t  y 


(7) 


where  shear  correction  coefficients  k and  k have  now  been  introduced 
in  a manner  similar  to  that  of  Timoshenko  [12]  and  Mindlin  and 
Deresiewicz  [17]. 

The  procedure  involved  in  deriving  the  theory  will  be  to  manipu- 
late the  left  sides  of  equations  (6)  until  they  are  of  fhe  form  of  the 
left  sides  of  equations  (7).  Thus,  taking  the  first  time  derivatives 
of  (6)  and  multiplying  by  the  viscoelastic  constant  C results  in  the 
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S.  P.  TIMOSHENKO  (1922)  Philosophical  Magazine,  Ser  6,  Vol  43, 
125-131.  On  the  Transverse  Vibrations  or  Bars  of  Uniform  Cross- 
Section. 

R.  D.  MINDLIN  and  H.  DERESIEWICZ  (1953)  Columbia  University  Tech- 
nical Report  No.  10.  Timoshenko's  Shear  Coefficient  fo**  Flexural 
Vibrations  of  Beams . 
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equatfans 


C ^y  *"  W 


C “y,y  + C °yz,Z  = pC  V 


(3) 


which  when  added  to  their  counterparts  in  equation  (6)  become 


°yz.y  + c V.y  = pM  + pc  ” 


°y,y  + C °y,y  + °yz,z  + C <W  = pV  + pC  v 


(9) 


Multiplying  the  first  equation  of  (8)  by  to  aria  the  second  equation  by  v 
integrating  over  the  beam  volume  and  time,  and  finally  adding  the  final 
answers  results  in  the  equation 


/ /V 

A o o 


(oyz,y  + c °yz,yW+  (V  + c ‘y.y^ 


+ (ayz,z  + ^ ayz,z^v 


dA  dy  dt 


■ i,  C C p 


V v + C *v  v + W w + C w*  w) 


dA  dy  dt 


(10) 


After  several  integrations  by  parts,  equation  (10)  may  be 

—i  % 


expressed  as 
,t 


J J 

A 0 


(o  + C a )w  + (a  + C 6 )v 
yz  yz  y y 


dA  dt 


0 


(°v,  + 0 aV7)v 


dA  dy  dt 


* s.  I-  /*  A 


masna* 


I I'j" 

A 0 0 


(ov  + C 6V7)  (—  + — ) 

v yz  uyz'  '.gy  3z ^ 

/ v 3v 

+(oy  + c V 5P 


dA  dy  dt 


= / 

A 0 0 


(w  + C w)  w + (v  + C v )v 


dA  dy  dt 


it  is  immediately  clear  that 


/ A A - 

A ° ° dz 


(cyz  + c V* 


dA  dy  dt  = 0 


since  both  beam  surfaces  are  stress  free  and  that 


OD 


(12) 
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j 

$ 
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/ / 

A 0 


(°yz  + C V)  * + (oy  + C °y)v 


dA  dt  = 0 
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since  the  boundary  terms  will,  be  satisfied  at  the  beam  ends.  ' Applying 
equations  (5)  and  (7)  to  equation  (11)  and  taking  into  account  equa- 
tions (12)  and  (13)  results  in 


J A A 

JQ  J0 


(2kGeyz  + 2k*G*eyz)(2eyz) 

•fa  * t 

+(ECy  + E Cy)ey 


dA  dy  dt 


/ ft  p[(W  + C w )w  + (V  + C v )v]  dA  dy  dt  . (14) 

A 


But,  from  the  chain  rule  of  partial  differentiation  it  is  clear  that 

8 r 2n  . . 

^Le  J = ee  + ec 


or  that 


(15) 


Similarly,  the  fact  that  an  indefinite  integral  can  be  defined  as  a 
definite  integral  with  a variable  upper  limit 

/g(t)dt  = /atg(t)dt  + const.  (17) 
immediately  results,  afte"  taking  a time  derivative  of  both  sides,  in 
the  equation 


~ /t  g(t)dt  = g(t) 
dt  a 

.2 

which  for  g(t)  - c results  in  the  relationship 

2 d t 
e = - / (e'r dt 

dt  o 


(18) 


(19) 


A direct  application  of  relations  (16)  and  (19)  to  equation  (14)  with 
an  introduction  of  equations  (4)  and  (5)  results  in  the  equation 


/*■  A d 


'0  '0 


3t 


1 AkG ( — - 4))2  + AkV/t(-  -i)2dT 

2 ay  ;o  ay 


AE(z^)2(^,  AE*/*(z^)2(||)2  dz 


* E(!l>2  ♦ Eh/hit,2  dt 


dy1 


dy  dt 


+ [£  ft  P _d 

* r\  Hi* 


AW2  - 2AC/0tw2dx  + A(z^)2i2 


,,  i + • 2AC/t(z.k)2rdT 

'o  o 2 at  I * ;o  a 


- 2IC  /Vdx 
Jo 


dy  dt  = o 


(20) 
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after  an  integration  over  the  beam  area  where 


A = bd 


with  b being  the  beam  width  and  d being  the  beam  thickness. 

Following  Anderson  [18],  a conservation  law  is  sought  in  the 
existence  of  a quantity  H such  that 


H = constant, 

(22) 

such  that 

obviously 

s ■ 0 

(23) 

where 

H = T + U + V 

(24) 

with  the  quantities  T,  U,  and  V being  called  the  kinetic  energy,  the 
potential  energy,  and  the  dissipation  energy.  From  a comparison  of 
equations  (20),  (23),  ar>  (24)  it  is  clear  that  the  various  energies 
may  be  defined  as 

T ■ £ £ T*  dy  dt 

li  ■ £ £ U*  dy  dt 

V = J*  J*  V*  dy  dt  (25) 

J\,Q 


G.  L.  ANDERSON  (1975)  Journal  of  Sound  and  Vibration,  39(1), 

55-76.  Stability  of  a Rotating  Cantilever’  Subjected  to  Dissipative, 
Aerodynamic,  and  Transverse  Follower  Forces. 
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and  from  equation  (20)  it  is  clear  that  these  energies  are 


T*  P rn  • 2 ^ V 2 * 2 T*2~, 

T = - [Aw  + A(za)  t|;  + 4 J 


u*  = I AkG(^  -*)2  <■  -(zk)2(^)2  + y<^)2 
2 Sy  2 a 3y  2 


*kV(g  4)2  ♦ e*k|)2 

+ AE*(zk)2(^)2  -pAC  «2 
a ay 

- pAC  (z£)2tf2-pIC  $2 


THE  LAMINATED  BEAM  THEORY 


The  laminated  beam,  Figure  1,  is  composed  of  a number  of  alterna- 
ting plane,  parallel  layers  of  two  homogeneous,  isotropic  viscoelastic 
materials  which  are  respectively  termed  the  reinforcing  layer  and  the 
matrix  layer.  The  reinforcing  layer  is  the  stiffer  of  the  two  layer 
combination  and  is  indicated  by  the  subscript  "1"  while  the  softer 
matrix  layer  is  indicated  by  the  subscript  "2".  The  elastic  constants, 
the  viscoelastic  constants,  the  layer  density,  and  the  thickness  for 
the  reinforcing  and  matrix  layers  respectively  are  E^  G-j , E^,  G-j , 

C-| , p i > d-j  and  E2>  G2,  E2 » ^2*  ^2’  ^2’  ^2’ 

The  basic  variables  involved  are  w,  the  transverse  deflection;  >J >,  the 
gross  rotation;  cj^ , the  rotation  of  the  stiff  layer;  and  cj^,  the  rota- 
tion of  the  soft  layer.  The  midplane  positions  for  the  kth  pair  of 

k k 

neighboring  reinforcing  and  matrix  layers  are  y-|  and  y2  respectively 


EE 


ii  in  !■  ii  i 


\ 


a?,  indicated  in  Figure  1,  with  the  layer  midplones  taken  perpendicular 
to  the  z-axis.  The  width  of  the  beam  is  b and  the  total  or  gross  thick- 
ness is  h. 

From  equation  (26),  the  kinetic,  potential,  and  dissipative  ener- 
gies in  tne  individual  layers  are 


Tak  * ♦ Aa(z*fy2  + laia2] 

C ■ \ Waify  V 

4.  !aEa/ 7k\2/3jtx2  + ^/^2 
2 [7*'  l3y;  2 9y 


V*k  = 
a Jn 


A kV(S  -i  )2  ♦ £*I  (~V 
a a a 3y  a a a Dy  ‘ 

+ A Et(zk)2(|^)2  -p  A C w^ 
a «'  a'  v3y'  a a u 

- p A C (zk)2if2  -p  I C $2 

Ot  Ct'  Ct/  r 1 n»  rvT/v 


a a a a 


dx, 


(2/) 


where  a r-  1,2  respectively  gives  the  reinforcing  and  matrix  layer 
energies. 

Now,  the  three  energies  are  summed  over  the  n layer  pairs  to 
determine  the  total  energies  for  the  composite  beam 
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* k=n  *1  *b 
T = l (Tk  + T*k) 
k=l  1 


u*  = kf  (u;k  ♦ ujk) 

k=l 


. k-n  ,A.,  *iy 

v - l (V,k  + v k) 
k=l  1 L 


m 


It  is  now  convenient  to  convert  the  discrete  system  (28)  to  a conti n-  t 

uous  system  by  utilization  of  a smoothing  operation,  that  is  to  replace 
the  summations  in  (28)  by  weighted  integrations  over  the  tnickness 
variable  z. 

The  result  of  the  smoothing  operation  is  the  energies 


* h/2  i 

T « [ 1 

-h/2  (dl+d2) 

* h/f2  1 

-h/2  (dl +d2^ 

. hf  1 

-h/2  .(d-t+d2) 


(T*  + Tp  dz 
(U*  + Up  dz 
(V}  + V*)  dz 


(29) 


where  after  smoothing 

z « z!f  ■ 7^  . (30) 

Carrying  out  the  integrations  in  (29)  in  terms,  or  (27)  and  taking 
into  account  (30)  results  in  the  energies 
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T‘  - jtp/,  + P2A2)(d^d7)  *2  ■''  2^(plAl  + P2A2)(?[Td7)  *■ 

1 h • 2 1 h • 2 

+ 2 Pl!l  (d^  'l’l  + 2 P2!2  (d^Tj  *2 

U‘  ‘ \ Wl  (^)(ry  V + l WZ  <d£i><®  V 


1 , „ r , a r \ h3  / , 1 P T h ,^1  N 

24  A1 E1  2 ElI]  (d1+d2)  3y  } 


3+i.2 


* 1 hh^ 


»*■/* 

0 


. ★ h / 3w  ’ >2  *p*  h / 3w  ',  \ 

A]klGl  Uh^diliy  ’ ^ A2k2G2(d1+d2)(3y 


1 * 

T2(A1E1  ' ”2-2'(d1+d2)'ay 


h3  /3iK2  . -*»  h ,3())1  \ 2 


+ ±&X  + *X)t7T7-A%)  + El"Il7Z^V) 


(d1+d2) 


* 

+ Eol 


l.  34>„  o 

l2l2TdT^d7T(3y“)  " (P1A1C1  ' P2A2C2)('d7^ 


h ..2 

w 


dr. 

(31) 


“ 12{P1A1C1  + p2A2C2)Tdj^  ^ “ °'1lCl((J1+d2)  $1 


h 2 
p9a9c97— - — - r 

M2  2 2(d1+d2)  '2 


At  this  point,  continuity  of  displacement  at  the  interface 
the  kth  pair  of  layers  must  be  considered.  Applying  equation  (l;  to 
a multilayer  beam  results  in  the  equation 

va(y,z,t)  = -z£  (y,t)  -z<j>a(y,t)  (32) 
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and  with  the  aid  of  Figure  2 it  Is  clear  that 

i<  . V 

V,  ■ -2,*  + -j*, 

k ^2 

V2  = ~Y  " 7*2  (33) 

at  the  interface  between  layers  1 and  2.  It  1:  also  clear  from 

i 

Figure  2 that 

zjj  ■ z*  - ^ + d2)  (34) 

and  that  equations  (33)  describe  the  same  interface  such  that 

v1  = v2  . (35) 

From  equations  (35)  applied  to  equations  (33)  it  is  clear  that  the 
continuity  condition  is 

iJj  = n<j>^  + 0-n)<f>2  (36) 

where 

di  do 

n ■ rt— h . (l-n)  = ?T~-t  . (37) 

(d*j+d2)  (d-j  +d2 ) 

Following  Sun  [8],  the  variable  <j>2  is  eliminated  such  that 


♦ 


2 


(l-n) 


(38) 


where  for  convenience  the  notation  <{>  = has  been  introduced. 

Expression  (38)  is  directly  substituted  into  equations  (31)  and 
the  dimensionless  variable 

C = h/(d1  + d2)  (39) 

8c.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  947-954. 
Microstructure  Theory  for  a Composite  Beam. 
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is  introduced  to  yield  the  energy  expressions 

T*  = ^(p-|A-]  + P2A2^  ^ + ^b^npl  + ^_ri^p2^ 


1».  T , U T r <t'  P-~4)2 


+ ;5piIi*  + 


2’  2 ?Ml-n)  (1-n) 


U*  . -*)2  * 1W2(|  - ^ + ^ ,)2 


+ irV,Ml-n)E2)(g)2  + lE1i1(»)z 


+ U ' (A)  f/ 


v*  - J* 

0 


rfl  i*r>*/8w  *\2  , r,  ,*r*/9W  ^ 2 

^A1  kl ^1  ( ay  " ^ rjA2k2G2^9y  ’ (l-n)  (l-nr 


+ IbtnE*  + (l-n)E2](^)2  + SEi>i<5$)2 


+ CE2l2(rr7o  55  " (A)  5yl2  '5(piAici  + P2A2C2)S 


- IjjtnPiC-,  + (1  -n)p2C2^2  " ^PlAA  $ 


dx 

(40) 


where 


lu  = 


bir 


b " 12“  . 


(41) 
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Now,  all  the  squares  of  the  various  sums  in  equations  (40)  are  expanded 
out  to  yield  the  final  forms  of  the  energy  expressions  as 


T*  - |Sa4w2  + jF,agj»2  + |sa13i2  -5a,,,# 


lr  /3Wv2  ,3w  , ‘\c  ,2  9w 

U = 2^al  Oy  "?a3*3y  2^a12^  '^a2l^y 


+ ^a6*2  -5a8#  + |cas(|^)2  * Ica,,^)2 


•k  .t 

v = jz 
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34/  34> 

K^Vy  ~y 


,3Wv2 


3w 


.2 


3w 


»/i2 


^1^)  -2f.b3*5-  +50, 2*  -25b2fe  + Cb6^ 


■ 3'K2 


,3^2 


3(Jj  3 if; 


■2*8**  + «b5(i>  + «bU<5p  dr  (42) 


f>4«2  _^bg>^  “ j ^ + 2£b,J| 


•{2 
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where  the  constants  a^  are 
a-j  = A i k -j  G i + 
a 2 = A2  ^2^2^ ( 1 -n) 
a^  = A^k^  -nA2k2G2/(l-n) 

a4  = d-jA-j  *!•  P2A2 

Ik  E?^2 

aR  = T-CnE-i  + (l-n)H2]  + 5 

551  2 (l-n)2 

j6  = A2k2Sz/(l-ri)2  = a^d-n)  i| 
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2 A2k2G2  = nafi 


p2!2 


aq  = —[npi  + (l-n)p?]  + 0 

9 « 1 2 (1-n)2 

n 
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all  " El'l+  ^2  E2!2 
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(1-n)' 
-.2 
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(1-n) 


p9 1/ 


2 M2A2 


(43) 


which  correspond  to  the  elastic  constants  given  by  Sun  [8]  for  elastic 


laminated  beams  and  where  the  constants  b^  are 


b1  - A^G*  + A2k*G* 


b2  = A2k*G*/(l-n) 


b3  = AlklGl  ' (1^)  A2k2G2 

^4  = Pl^iC-j  + P2^2^2 
Ii 


ih  . Eolo 

b5  - ~[nE*  + (i-n)E2]  + ^-L 


(1-nV 


b6  = A2k2G*/(l-n)2  = b2/(l-n) 


i_  _ n ,-*T 
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C.  T.  SUN  (1971)  Journal  of  Applied  Mechanics,  38,  947-954 
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b8  = A2k2G2  = nb6 

bg  * ^LnP1C1  + (l-n)p2C2]  + p2I2C2 


HO  " 2 P2I2C2 

(1-nr 


* 2 
* n r* 


b = E-]  I -j  + -y—  1 2 

11  11  (1-n)2 


* * 


bio  = A.k.G,  + — A_k„G 


* * 


12  "m 


(1-n) 


2 2 2 2 


bl 3 = P1I1C1  + P2I2C2  (44^ 

which  correspond  to  the  viscoelastic  contribution  of  the  current  anal- 
ysis of  viscoelastic  laminated,  beams.  It  should  be  noted  that  the 
author  [19]  has  evaluated  viscoelastic  shear  correction  constants  in 
another  paper  and  based  on  this  evaluation  it  is  clear  that  k-|  = 
k2  = k*  = kj  = it2/!  2. 

Now,  from  equations  (22-25)  in  conjunction  with  equations  (42)  it 
is  easy  to  form  energy  principle  (23),  that  is  dH/dt=0,  which  upon 
various  integrations  by  parts  and  a gathering  of  common  factors  of 
w.ijj,  and  $ results  in  equation  (23)  becoming 
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dt  = 0 . ' (45) 
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The  viscoelastic  equations  of  motion  and  boundary  conditions  for 
laminated  beams  are  now  obtained  by  applying  the  first  'lemma  of  the 
calculus  of  variations  to  equation  (45).  Thus,  the  throe  equations  of 
motion  are 

32w  S'!1  3<j>  , , 32w  , 3<i>  I 3d)  ..  , ... 

al^  “ a23y  _ a33y  + bl^  " b23y  " b33y  = a4w  + b4w 


ay 

32ci. 


3w  32t|j 


a2i  + *5^  - V - a7^f  + *8*  + % + b5~2  - V 


3y 


,2 


by— t + bgii  3 agip  - a-]Q$  + bg«j>  - b-|0^ 


ay 

a - . a ^ + a i|i  + a ^ - a <j)  + b — - b^’-<bij> 
33y  73y2  8*  ll3y2  1 - 2 ay  7^  8^ 

+ b1irf  ‘ b12*“'al0if  + al3;f-b10*  + b13* 

ay 

and  the  corresponding  boundary  conditions  are 
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y = 0,A 
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Ip  = 0 on  y = 0,£. 
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or 


(()  - 0 on  y = q,sl 


(47-c) 


SUMMARY 

An  energy  principle  has  been  formulated  for  viscoelastic  Timoshenko 
beams  according  to  the  standard  linear  model  with  the  stipulation,  and 
hence  additional  terms,  that  the  energy  principle  be  utilized  in  build- 
ing a viscoelastic  laminated  beam  theory.  The  Timoshenko  model  con- 
sidered has  accounted  for  both  viscoelastic  extensional  and  visco- 
elastic shear  strains.  To  later  incorporate  the  single  layer  energy 
principle  into  the  development  of  a laminated  Deam  theory,  a term  which 
accounts  for  the  beam's  gross  rotation  was  included  in  the  single  layer 
development. 

Using  the  single  layer  energies  developed,  a viscoelastic  lam- 
inated beam  theory  composed  of  a number  of  alternating,  plane,  parallel 
layers  of  two  homogenious,  isotropic  viscoelastic  materials,  termed 
the  reinforcing  layer  and  the  matrix  layer,  was  derived.  In  deriving 
the  theory,  the  individual  layer  kinetic,  potential,  and  dissipative 
energies  were  summed  over  n layer  pairs  to  obtain  the  total  energy  of 
the  compos i te  beam;  these  results  are  converted  to  a continuous 
system  by  utilization  of  a smoothing  operation  or  weighted  integration. 
The  number  of  independent  variables  in  the  total  composite  beam  energies 
is  reduced  from  four  to  three  thru  the  introduction  of  a condition  for 
continuity  at  layer  interfaces;  A direct  application  of  the  energy 
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principle  developed  to  the  composite  beam  energies  results  in  a set  of 
three  equations  of  motion  and  their  corresponding  boundary  conditions 
for  viscoelastic,  laminated  composite  beams.  A future  report  will  dis- 
cuss a number  of  meaningful  results  of  both  wave  propagation  and  vibra- 
tion analysis  utilizing  the  equations  of  motion  developed  in  this  report. 
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APPENDIX 


The  present  objective  is  to  derive  a set. of  constitutive  relations 
which  can  be  utilized  in  conjunction  with  the  basic  equations  for  a 
Timoshenko  beam.  While  constitutive  equations  may  be  formulated  in 
either  integral  or  differential  form,  preliminary  work  in  the  direction 
of  formulation  of  a viscoelastic  beam  theory  for  laminated  composite 
materials  indicates  that  the  differential  form  of  constitutive  rela- 
tions will  be  most  useful.  The  differential  constitutive  relations 
will  be  utilized  in  the  present  development. 

The  general  form  of  the  differential  constitutive  equations  is 
adapted  from  Fung  [20]  where  the  stress-strain  relations  are  of  the 
form 

" V0)eij 

P2(D)okk  ■ Q2(l))ekk  (fl-D 

where  P^(D)  and  (D)  are  given  by 

k=n-i  . 

I’,(D)  = T akDk 
k=o 

k-n,  k 

P2(D)  = l 2CRD 
k-:o 


20 

Y.  C.  FUNG  (1965)  Prentice-Hall,  Inc.  "Foundations  of  Solid 
Mechanics" . 
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and  operators  (A-2)  in  light  of  (A-6)  reduce  to 
P-]  (D)  = aQ  + a-j  D 

Q-j  (D)  = bQ  + b-,D 

P2(0)  = C0  + C1  0 

Q2(D)  - d0  + d,D  . (A-7 ) 

As  will  be  subsequently  seen,  the  only  non-zero  stresses  and 

strains  for  a Timoshenko  beam  with  its  y-axis  along  the  length  and  its 

z-axis  through  the  thickness  are  a and  a & e and  e , Thus,  from 

y y^  y y ^ 

equation  (A-4)  the  non-zero  stress  and  strain  deviators  are 
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(A-8) 


Now,  a direct  substitution  of  equations  (A-2),  (A-7),  and  (A-8)  into 
equation  (A-l ) results  in 

[1  + (a-,/a0)D]ay7  = [(b0/aQ)  + (b1/aQ)D]eyz 
[1  + (a1/a0)D]ay  = [(bQ/a0)  + (b1/aQ)D]ey 

[1  + (C1/C0)D]ay  = [(d0/CQ)  + (d1/CQ)D]ey  . (A-9) 

There  are  thus  two  equations  for  stress-strain  in  the  y-  coordinate 


Dl°y  = °2ey 
°3uy  = D4ey 


(A- 10) 
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where 


O')  = l-(a-|/a0)D 

°2  = (bo/ao}  + (b1/ao)D 
°3  = 1 + (C1/Cq)D 

°4  = (d0/c0)  + (dl/c0)°  * (A-ll) 

and  they  must  be  combined  to  form  a single  constitutive  equation 

20^30  « (D2D3  + D] D4)ey  . (A-12) 

Now,  from  both  the  right  and  left  sides  of  equation  (A-12)  it  is  clear 
that  the  constitutive  equation  is  of  the  form 

(1  + aD  + bD 2)oy  = (1  + C D + dD2)ey  , (A-13) 

but  it  would  now  be  desireable  to  have  the  form  of  the  standard  linear 
model  as  in  equation  (A-5),  if.  possible.  This  can  be  achieved  if  the 
restriction  is  now  made  that 

D1  = D3  = 1 + (a1/aQ)D  (A-U) 

such  that  equation  (A-12)  now  becomes 

[1  + (a,/a0)D]oy  - l[(b0/a0  ♦ d0/C0)  + (b,/a0  + <yC0)D]Ey  .(A-15) 
As  a final  step,  define  the  constants 

C = a-j  / a0 

E " 2(bo/ao  + ^0^0) 

e‘  - 5(n/a0  + d,/C0) 

2G  > b0/a0 

2G*  = b./a 
1 0 


(A— 1 6 ) 


with  the  final  form  of  the  constitutive  equation  thus  being 
(1  + C jpj:)oyZ  = (2G  + 2G  ^)ey2 
(1  ♦ C ±)ay  = (E  ♦ E*j|)ey  . 


( A- 1 7 ) 
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